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ABSTRACT 

A solution of the large discrepancy existing between inclusive and exclusive measurements of 
the 8 Li + 4 He — > n B + n reaction cross section at E cm < 3 MeV is evaluated. This problem 
has profound astrophysical relevance for this reaction is of great interest in Big-Bang and r- 
process nucleosynthesis. By means of a novel technique, a comprehensive study of all existing 
8 Li + 4 Hc — ► n B + n cross section data is carried out, setting up a consistent picture in which 
all the inclusive measurements provide the reliable value of the cross section. New unambiguous 
signatures of the strong branch pattern non-uniformities, near the threshold of higher n B excited 
levels, are presented and their possible origin, in terms of the cluster structure of the involved 
excited states of n B and 12 B nuclei, is discussed. 

Subject headings: early universe — nuclear reactions, nucleosynthesis, abundances — supernovae: general 



Recently, the total disagreement between in- 
clusive measurements of the 8 Li+ 4 He— > n B+n 
reaction cross section and exclusive ones at kinetic 
energy in the centre-of- mass system E cm < 2 MeV 



has been pointed out (|La Cognata et al. l2008h . 



The case is illustrated i n FigED The comp l emen 



tary neutron inclusi ve (ILa Cognata et al.l 12 008) 



and 11 B inclusive ( Bovd et al. 1992t Gu et al 



1995) measurements give comparable values of 



1 Supported in part by Deutsche Forschungsgemein- 
schaft. 



the reacti on cross section . Inst ead, the exclusive 
approach l lshivama et al. 2006), where 11 B-n co- 
incidences were measured, gives cross-section val- 
ues which, with the exception of a narrow region 
around E cm ~ 1.7 MeV, are generally smaller by 
a factor < 4 with respect to the inclusive ones. 

The largest discrepancy among the data in 
Figjl] is observed right at energies of astrophys- 
ical interest E cm = 1 MeV. The measurements 
at such energies explore the Gamow window of 
the 8 Li+ 4 He^ n B+n reaction at temperatures 
T = (1 -T- 3) X 10 9 K, of great interest in primor- 
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dial as well as in other relevant nucleosynthesis 
sites, in particular corc-collapsc supcrnovae and 
neutron-star mergers. Within the frame of the in- 
homogeneous Big Bang model, still rep resenting a 



viabl e possibility for the early un iv erse (ILara et al 



2006 



Malanev fc Fowler 



1990t iRauscher et al.l 120071 ) , this reaction could 



1988; iKaiino fc Bovd 



have allowed to overcome the ^4=8 mass gap, thus 
providing a possible explanation for the experi- 
mental observation of a non-negligible abundance 
of he avy elements in the old est astrophysical ob- 
jects ( Matsuura et al. ( 2005f ) and Refs. therein). 
The magnitude of the cross section is the key infor- 
mation to check the reliability of these predictions. 
The 8 Li+ 4 He^ 11 B+n reaction also plays an im- 
portant role i n the context of the r-process nu- 
cleos ynthesis ( Sasaaui et al. 20061 : Terasawa et al. 
l200ll ). Currently, the most popular scenario is 
neutrino-driven winds from Type II SNe. Any- 
way, possibility remains that it could be asso- 
ciated with neutron-star mergers or gamma-ray 
bursts, in which the required neutron-rich con- 
ditions can also be realized. Therefore, it is of 
critical importance to constrain the parameter 
space for t he r-process to restri ct possible envi- 



Sasa qui et al. (2006) found that the 



ronments 

8 Li+ 4 He— > 11 B+n reaction leads to a more effi- 
cient production of seed nuclei, so that a larger 
neutron/seed ratio is required for a successful r- 
process. This, in turn, allows to constrain the 
entropy per baryon and the astrophysical site for 
production of r-process nuclei. 

In this Letter, we conclude that inclusive mea- 
surements provide the most reliable estimate of 
the cross section. 

For the 8 Li+ 4 He— > 11 B+n reaction involves the 
unstable 8 Li nucleus (T1/2 = 840 ms), the scope 
of measurements with present-day facilities is the 
cross section a(E cm ) summed over all 11 B+nj 
branches, 



o-(Ecm) 



E 

i=0 



N 



del 



^tar^proj 



(i) 



imax denoting the energetically open highest 
branch at a given E cm (i = being the branch 
leading to n B ground state). The involved re- 
action bra nches, identified by th e 11 B level se- 
quence 
TabED 



Ishivama et al.l ( 20061 ). are listed 



might spoil cross section measurements, the num- 
ber of detected particles N^ et , of impinging pro- 
jectiles N 1 "" -?, of target nuclei per unit surface 
N tar , and the detection efficiencies £j contribute 
to the overall uncertainty. Background might illu- 
sorily enhance t he reaction yields et of inclusive 
measurements (Boyd et al.l Il992t iGu et al. 1995: 



La Cognata et al.l 2008). However, they were per- 



formed by measuring completely different ejec- 
tiles ( n B or neutron), thus different sources of 
background are expected. Therefore, the agree- 
ment between them make us confident that the 
possible background contribution stays below the 
uncertainty ranges in Fig[T] We also underscore 
that each of inclusive dBovd et al.lll992l:lGu et al ' 



1995 ) and exclusive ( Ishivama et al. 20061 ) data 



sets has been measured in a single irradiation, 
thus any error in their evaluation attributable 
to N tar and/or W roj only could lead to a rigid 
vertical displacement of one excitation function 
with respect to the others, e xcluded by Fig[T ] 
Indeed, concor dant inclusive ( Bqvd etabl 199 



Gu et all Il995h and exclusive (jlshivama et al 



20061 ) cross section measurements at E c . m . ~ 1.7 



McV (sec FigQ} rule out significant errors on both 
N tar and N pr ° : ' . Even supposing the concurrence 
of equally oriented systematic uncertainties on 
Nf, N* ar and W roj it is not possible to ex- 
plain the magnitude of the discrepancy and its 
dependence on E cm . The most likely candidate 
source of error is then represented by the neu- 
tron d etection efficiency in t he exclusive measure- 
^ 20061) 



ment ( Ishivama et al 



Indeed, in compar- 



ison to B inclusive measurements ( Bovd et al 



199i IGu et al.lfl995h . the sensitivity to reaction 
n exclusive measurements 



events in the B 
(|lshivama et al. 2006) was governed by the neu- 
tron counter. In comparison to neutro n inclusive 
measurements ( La Cognata et aLll2008l ) asubstan- 
ti ally different type of n eutron detector was used 



lly dmcrc nt type or n eutron detector was used 
Ishivama et al.l ( 20061 ). As s hown in FigfJl the 



thcrmalization counter used in La Cognata et al 



Examining possible error sources that 



(|2008l ) is a zero-energy-threshold detector, its de- 
tection efficiency staying at a significant level 
down to thermal e nergies. Instea d , the plastic 
scintillator array in llshivama et~aH (|2006l ) shows 
a steep drop in detection efficiency with decreasing 
neutron energy E n below 2 MeV with a seeming 
cut at E n = 0.5 MeV. The occurrence of such 
a detection threshold plays a critical role in the 
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measurement of the total cross section according 
to EqQ] With this respect, we consider the calcu- 
lated kinematical diagram for the n B final states 
of the 8 Li+ 4 Hc^ n B+?i reaction at E cm = 0.75 
MeV, shown in Fig[3J In this plot, with the cut 
at 0.5 MeV in the neutron energy spectrum, the 
i = 6 reaction branch is experimentally unacces- 
sible. In such a situation the exclusive exper- 
iment ( Ishivama et al. could not provide 
the wanted cross sections summed over all 11 B 
final sta tes. Moreover, w i th the i nformation avail 



able i n llshivama et al. (l2006h : lHashimoto et al 



(120061) we have performed Monte Carlo simu- 
lations by imp le ment in g the detector se t- up in 



Ishivama et "ail (|2006l ); [Hashimoto et all (120061) 



into a GEANT code as described in Celano et al 



(|I9970 . Excellent agreement with the bulk of the 
experimental efficiency curve is obtained by the 
dashed histogram in FigEl This shows a steep 
drop below 2 MeV, very similar to a sharp de- 
tection cut-off, characterized by the effective half- 
drop energy threshold E^ ut = 1 MeV. With refer- 
ence to FigEl a cut-off of about 1 MeV would make 
unaccessible the branches i = 4, 5 besides the i = 6 
reaction branch. Quantitatively, the more favored 
the associated branching ratios the larger is the 
missing cross section in the exclusive experiment. 
It should be noted that the case E cm = 0.75 MeV 
illustrated in Fig[3]is in the hearth of the Gamow 
window for Big Bang nucleosynthesis. 

The laboratory reaction neutron kinematics is 
summarized in Fig [4^ as a function of E cm . For 
each reaction branch, at a fixed E cm , the labora- 
tory neutron energy ranges between the minimum 
Ef" n (dashed) and the maximum i^ max (solid) 
curves. Assuming E™* = 1 MeV, we have calcu- 
lated the E cm values t^ ut at which £™ ax crosses 
such E™* level (dashed black line in Fig|4jt). 
These are given in TabQ] together with the cor- 
responding set of reaction threshold energies U 
at which a branch starts to be potentially ac- 
tive, fixed by the corresponding reaction Q- values. 
Each ti — i?"' couple singles out a E cm kinematic 
region where the observation of the corresponding 
i— th branch is completely missed, because the cor- 
responding laboratory neutron energies stay below 
the experimental threshold. These ti -f- tf ut inter- 
vals are emphasized by shaded bands in Fig|4^. 
Inside each E cm shaded interval, no efficiency cor- 
rection can by any means be performed for the 



corresponding completely missed branch. 

Whenever E cm falls outside the shaded region 
of an open branch this becomes partially observ- 
able and its missed portion is recoverable, only 
provided that the observed portion is corrected 
by the appropriate value of its observability fac- 
tor P[ ut . For each reaction branch, this is de- 
fined as the observable portion of the laboratory 
neutron energy distribution -j^- characterized by 

Consequently, in the exclusive measurements 
( Ishivama et al.ll2006l ). the excitation function, i.e. 
the cross section versus E cm , can be corrected at 
some E cm -values, but remains necessarily uncor- 
rected at some other E cm values, right inside each 
shaded interval, where the P[ ut of the correspond- 
ing i = 6-^9 reaction branch equals zero. 

Conversely, the zero-energy detection thresh- 
old inclusive measurements (|Boyd et al. 119921: 
Gu et all [l99l lLa Cognata et al.l 120081 ) do 



mea- 



sure the wanted cross section a in the whole in- 
vestigated range of E cm . Accordingly, we write 



l— 2 ma> 

E 

i=0 



(2) 



where the factors fi are the unknown branching 
ratios as functions of E cm . In EqJ^l (Tind is the 
weighted line ar interpolation of all in c lusive cross 
section data ( La Cognata et al.l 120081 ; iBovd et af 



I992t iGu et al.lll995l ) within 0T MeV bins, repre- 
sented by the curve in Fig[T] 

On the other hand, in the case of t he exclusive 



20061 ) the detec- 



measurcmcnts (jlshivama et al. 
tion efficiencies ti must incorporate the observ- 
ability factor P[ ut of the corresponding branch 
versus E cm . Assuming isotropic neutron emis- 
sion in the centr e-of-mass system, foll owing the 
measurements in llshivama et al.l (|2006l) . we have 
calculated the Pf ut values for different E 1 ™*. 
The results for E™* = 1 MeV arc displayed in 
FigEb- We remark that with the higher threshold 
-E™' ~ 1 MeV suggested by Figj2] the portions 
of missed even t s in t he exclusive experiment by 
(l2006lh 



Ishivama et al 



and hence the corrections, 
are much larger than those corresponding to the 
seeming cut of 0.5 MeV so that a P . X ri data have re- 
maine d essentially uncorrected in llshivama et al.l 
( 20061 ). Accordingly, the experimenta l cross sec- 
tion a exc i reported in llshivama et ah (2006) can 
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be written as 



Texcl = u- < P cut > 

t= tmax 

<Pcut>= J2 f*- P t Ut 



(3) 
(4) 



< Pcut > being the average of the observabil- 
ity factors P? u * weighted by the branching ratios. 
Eqs[3]|4] not only clarify the primary role played 
by the experimental threshold E™ 1 through the 
factors P[ ut (< 1), but also the fundamental con- 
current role played by branching ratios. Indeed, 
we underscore that if all Pf ut = 1, a exc i = a 
no matter the branching ratio pattern because 
EilJ m " fi = 1 an d, consequently, < P cut >= 1. 
Conversely, if at least one of the P(' ut is sig- 
nificantly smaller than one a exc i underestimates 
cr, the more favored the feeding of a (partially) 
missed branch, the larger the deviation of <J exc i 
from a. With reference to the case in Fig[3l this 
effect is more significant when reaction branches 
involving large 11 B excitation energies come into 
play, their observation being extremely sensitive 
to the neutron threshold energy . Clearly, 
the deviations of <r exc i from <Ji nc i are entirely de- 
scribed by < P cut >. Accordingly, from EqUH we 
have deduced the exper imental < P ru t > of the 
exclusive measurement ( Ishivama et al.l 120061 ) as 

< Pcut > = f 23 ^, the errors of the interpolated 
Pind being appropriately propagated. It is shown 
in FigHJ: versus E cm and confirms that the largest 
deviations between a exc i and ainci occur right in 
correspondence of the shaded bands, forming a 
marked saw-tooth-like behavior with two appar- 
ent falls corresponding to the opening of the i = 6 
and i = 8 branches. 

To describe such a rise-and-fall behavior, a sim- 
ple recurrence formula can be used, following im- 
mediately from Eq01 



<p c ( + } >=(i-/ im „j-<p c y >+fi 



■pf 



(5) 

where (— ) and (+) denote near-threshold E cm val- 
ues on the left and right side of ti mate , respectively 
The two terms on the right side govern the magni- 
tude of each fall and of the following rise, respec- 
tively. 

As a reference case, in FigHJ; we show the 
< Pcut > calculated according to Eq0] for the 
iWw* = {imax + 1) uniform branch pattern 



(thin solid line). The energetically open fi's are 
assumed constant in each U -r- ti+\ E cm interval 
and the P[ ut in Figib, evaluated for E^ ut = 1 
MeV, are used in the calculation. In this case, 
falls of < P cut > occur at the successive opening 
of each of the i — 4 -j- 9 reaction branches because 
the threshold E^ ut on the neutron energy makes 
each of them experimentally unaccessible inside 
the corresponding ti -f- f™' interval. Strikingly, not 
only the falls in the experimental < P cut > recall 
this type of discontinuity but, in addition, we can 
conclude that in the experiment the i = 6 and 
i = 8 branches strongly deviate from an uniform 
pattern. 

Because of the normalization to 1, when the 
i = 6 and i = 8 branches strongly add up, 
correspondingly strong falls should be observable 
by examining the population of all the lower-i 
active branches. Indeed, this is what we have 
found in the trend of the experimental /q = 
J° - values as a function of E cm , as it is clearly 
demonstrated in Fig[4ji. This has been deter- 
mined hCTe_s_tajting both from the exclusive do 
data ( Ishivama et al.l [20061 ) and from the ctq de- 



duced from the ll B(n , a 
(Il990h . 



Paradellis et al 



s Li inverse reaction in 
We also note that an 
apparent increase of fa in both data sets distin- 
guishes two regions, below and above E cm ~ 1.6 
MeV, where different branch pattern regimes come 
presumably into play. 

Accordingly, two separate fits of the experimen- 
tal < P cu t >, for the region below and above 
E cm ~ 1-6 MeV, has been performed to deter- 
mine the branching ratios fi, using the subrou- 
tine MINUIT. Following Fig[4ji, the branching ra- 
tios fo-Hmax are assumed constant inside each E cm 
interval in between two successive branch open- 
ings. Therefore, all /, can be treated as free fit- 
ting parameters only constrained by the normal- 
ization condition X)i=(T fi = ^- ^ n particular, 
as the factors Pq1\ = 1 (Figgb), EqU becomes 

< P cut >= J2to fi + E*=4 X fi ■ p t ut so that we 
have considered the i = -f- 3 as a single branch. 
Because the threshold energies £4 and £5 are so 
close (Tab[T]), also the i = 4 -=- 5 are treated as 
a single branch. The resulting branching ratio 
(and error) values are listed in TablTland the as- 
sociated curves are shown in FigdflJ. We also 



1 The sensitivity to each fi is the higher the steeper is the rise 
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remark that both start energy and shape of the 
two rises in the experimental < P cu t > are ex- 
tremely sensitive to By tests performed, 
this fact leads to the independent determination 
of E^ 1 ~ 1 MeV, in agreement with the analysis 
of efficiency data in FigfS] To cross check these 
results we evaluate the < ff ut > averaged over 
the 0.75 < E cm < 2.55 MeV rang e , and compare 
with those given in llshivama et all ( 20061 ) (Tab[T]). 
These < ff ut > arc linked with the true fi by 



/are °- • p r ut ■ f* ■ dE ° 



<ft 



>= 



fi ■ dE c 



T f 2 ' 55 cr • P c 
JO. 75 ■ 



(6) 



Inserting into Eq[6]the /, values established above 
results in the < ff ut > reported in Tab[T] the 
errors of a and of all the fi being propagated ac- 
cordingly. The good agreement of this comparison 
strongly supports the validity of EqJ5] and, there- 
fore, the reliability of <7j rac ;. 

The physical novelty of the present Letter is 
the apparent selective feeding of the highest ex- 
cited 11 B levels, in particular of the i = 6 and 
i = 8 branches leading to the n B levels at 7.29 
and 8.56 MeV, respectively, originating the dis- 
cussed non-uniformities. This peculiar trend can 
hardly be understood invoking selection rules only, 
as there does not seem to be anything unique in 
terms of quantum numbers about the i = 6 and 
i = 8 branches (see Table [T]). Rather it seems 
to signal that the nuclear structure of the initial 
12 B* and of the final 11 B* excited states plays the 
most important role in determining the charac- 
teristic non-uniformity of the observed branching 
ratios. Concerning 12 B, in the excitation energy 
region explored here (E* = E cm + 10.01 MeV), 
states with large a, t or 5 He spectroscopic factors 
have recently been emphasized in the 9 Be + 7 Li — > 
2a + 8 Li reaction. In particular, clear evidence 
exists for at least two states, at 10.9 and 11.6 



MeV , which show significant a widths (jSoic et al 
2003). The relative contributions of these 12 B 
states to the 8 Li(a ,n) 11 B reaction process are 
regulated by both their a and n partial widths. 
The 12 B state at 10.9 MeV clearly contributes, 
giving rise to a resonance at E cm ~ 0.9 MeV, as 
it is demonstrated by Oi nc i (FiglJ)- The contri- 



of the corresponding P? . This implies that in the E cm 
region, where the P,j™g slowly approach unity, the fitting 
procedure presumably adds small /4_6 values to 



bution of the 11.6-MeV 12 B state cannot be ex- 
cluded as a hump does appear right at E cm ~ 1.6 
MeV in the cross section (Ji nc i, indicating a small 
n-width and, therefore, a possibly complex clus- 
ter structure. Concerning n B, a well developed 
2a+t cluster structure of the i = 8 state at 8.56 
MeV has been established very recently, whereas 
the other lower- lying negative-parity n B levels 
are suc cessfully described by shell-model calcu- 



lations (jKawabata et al.l 120071 ). Right below 



the trend of the experimental < P cut > is consis- 
tent with the one obtained for the uniform branch 
pattern, which suggest an at most weakly non- 
uniform preferential feeding of the lowest energy 
i = -j- 3 reaction branches. This can be likely 
attributed to a mismatch between initial- (clus- 
ter) and final- (single particle) state nuclear struc- 
tures. Above ts, the enhanced relative feeding of 
the newly open i = 8 branch, signalling a large 
overlap between the initial and final state wave 
functions, indicates a 12 B cluster structure close to 
the 2a+t one of the 8.56 MeV level of the daughter 
11 B nucleus (plus neutron). 

The comparable reaction branch non-uniformity 
type established here for the i = 6 branch sug- 
gests close structures of both 11 B and 12 B involved 
states. The relevant issues pointed out here call 
for further investigations. 

In conclusion, in this Letter we have shown 
robust evidences that support the large values 



from inclusive measurements (ILa Cognata et al 
20081 iGu et~ai]|l995t iBovd et al.lll992f ) as repre 



senting at present the most reliable estimate of the 
8 Li+ 4 He— > 11 B+n reaction cross section at astro- 
physical energy The recommended value of this 
cross section (and the corresponding error Act), 
following the considerations developed above, is: 



a(E cm )=J2a k E^ m + J2 



k=0 



^{E cm -E 3 y+g*lA 



(7) 



Aa = J2c n El\ n (8) 



where a = -9.09 mb, a x = 97.2 mb/MeV, 
bi = 18.6 mbMcV 2 , b 2 = -514 mbMeV 2 , 
b 3 = 722 mbMeV 2 , g x = 0.509 MeV, g 2 = 
0.914 MeV, 33 = 1.05 MeV, E 1 = 0.986 MeV, 
E 2 = 1.84 MeV, E 3 = 1.87 MeV, c = 66.9 mb, 
ci = -18.7 mb/MeV, c 2 = 11.5 mb/MeV 2 . The 
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previous formula is a simple fitting of the cross- 
section weighted linear interpolation in FigJIJ 
whose accuracy is better than 5% in the whole 
energy range. 

The original approach we have developed here 
can have important applications in different fields 
and can be extended to become an effective ex- 
perimental method to extract spectroscopic in- 
formation otherwise inaccessible with present-day 
experimental facilities. For the nuclear physics 
case considered here we have determined for the 
first time significantly non-uniform branch pat- 
terns, which are interpreted as manifesting the 
exoti c cluster structure re cently dis covered in 



"B ( Kawabata et all 120071 ) and 12 B ( Soic et al 



2003) excited nuclei. More importantly, this Let 
ter strongly calls for revised calcu lations of the 
r-proc ess nucleosynthesis. Indeed, ISasaqui et al 



( 20061 ) concluded that the entropy per baryon in- 
creases by about a facto r of 2 from previous es- 
Sasaaui et a l. ( 2005), using the cross 



timates in 

section in Ishivama et al. ( 20061 ). According to 
the present revised cross-section value of the 
8 Li+ 4 He— > 11 B+n reaction, the consequent con- 
straint on models of the r-process astrophysical 
site might be significantly altered, with undoubt- 
edly interesting consequences for astrophysics. We 
will explore possible additional implications of this 
work in future studies. 
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Fig. 1.— 8 Li + 4 Hc 
section d ata versus E, 
• (iBovdetal 



L B + n reaction cross 



2008). 



▲ (jlshivama et al 



□ (ILa Cognata et al 



- 19921) . ■ (|Gu et alJIlOoli r 

2006|). The curve depicts the 



interpolated inclusive cross section <7i nc i- The 
Gamow energy region at T = 2 x 10 9 K is shown. 



7 




Fig. 2. — The dependence of the total detec- 
tion efficiency of the t hermalization counter in 
lLa Cognata et al. ( 2008T ) (solid his togram) and of 
the p la stic scintillator arrays in llshivama et al.l 
(|2006l ); iHashimoto et all (|2006f) (closed triangles) 
on the laboratory neutron energy. The dashed his- 
togram is the Montecarlo simula tion which better 



repro duces the efficiency data ([Hashimoto et al 
20061) . 
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10 




laboratory angle (deg) 



Fig. 3. — Kinematical plot showing the neutron 
energy-angle correlation for each branch of the 
8 Li + 4 He -> ll B + n reaction at E cm = 0.75 MeV. 
The minimum and maximum neutron energies al- 
lowed by kinematics in the whole angular range 
are given for the i = branch ( 11 B ground state). 
The red line represent the 0.5-MeV seeming cut 
on the neutron energy from FigO 
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(MeV) 



Fig. 4. — a) Kinematical loci of the minimum 
(dashed) and maximum (solid) laboratory neutron 
energies; b) the branch observability factors P£ ut 
for E™* =1 MeV. The curves of the branches 
i = 0^3, not shown above, are entirely above 
the E^ 1 = 1 MeV level and the corresponding 
pcut _ y independent of E cm ; (c) the experimen- 
tal < P cut > (symbols), its best fit (thick curves) 
and the uniform pattern case (thin curve); (d) 
the branching ratio fp deduced from groun d state 
cross sectio n data in Ishivama et al.l ( 20061 ) (open 
circles) and Paradellis et al.1 ( 199Clh (filled squares 
x2). The horizontal lines indicate the weighted 
mean fo values in the corresponding intervals. 



Table 1: n B final state, center-of-mass reaction energy threshold t, observability energy threshold t cut (assuming £yj"' = 1 MeV), true 
branching ratios / in the indicated E cm intervals and mean branching ratios < / >. Short lines indicate branching ratio values 
compatible with zero. Empty entries correspond to branches which are not physically open. In the fitting procedure, the i — ~ 3 branches 
are not separated and the summed /o-j-3 is given. The same applies to the i = 4, 5 branches. 
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